TWO-PHASE  LAMINAR  BOUNDARY  LAYER  FLOW  AROUND 

A  WEDGE 


Dirk  van  Dord 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


TWO- 

PHASE  LAMINAR  ] 
FLOW  AROUND 

BOUNDARY 
A  WEDGE 

LAYER 

by 

Dirk  van 

Dord 

June  1982 

Thes 

is  Ad 

visor 

James 

V. 

Sanders 

Approved  for  public  release;  distribution  unlimited 


-fi 


UNCTiA^S.'ST^Tm 


SCCUMITY   CUASSiriCATIOM   Of    THIS  m  tkOK.  (Wttat  Data  Bnl»t»d) 


REPORT  POCUMEHTATIOH  PAGE 


2.   OOVT    ACCKSSION   NO 


4      Ti  r\.C  fmt^  Subiltim) 


Two-Phase  Laminar  Boundary  Layer  Flow 
Around  a  Wedge 


7.      AuThOWiO 


Dirk   van   Dord 


i      *CMro«MINC  OAOi^NlZATtON    NAME   ANO   AOOACK 

Naval    Postgraduate    School 
iMonterey,    California      93940 


II     coNmoLLiNS  ornce  namc  amo  aoomcss 

Naval  Postgraduate  School 
Monterey,  California   93940 


14      MONlTOAlMS   AGENCY   MAMC   ft    AOOMCtSff/  tf<//ar«Mf  Ir^m  CaMlralXnf  OUIfi 


READ  tNSTR'jrT'ONS 
BEFORE  COMPLET:nc  FORM 


»     *eci»iSH  T's  c  A' *uoc  NijMSf 


5.      TY«»e    Of    «EOo«'    4    »£BioO   COVEPES 

Master's  Thesis 
June  19  82 


«.    ^mfomuina  o*g    mm 


0*T    NUMaCH 


•■    CONTRACT  o«  gbant  N».Mac«ri/ 


10.      »«OGaAM    ELEMCHT    BBOJECT       TASit 
AaCA   ft    «0«K   UNIT   NUMBCns 


la       HCPOXT    DATE 

June    198  2 


I)       NUMSCM  OF    "AGES 

62 


'»•     SeC'jBITY   CLASS,   (ot   (Ala  rtport) 

Unclassified 


'»•■      DECLASSIFICATION/  OOWNGRAOING 
SCHEDULE 


IS.     OlSTWiauTlON   STATEuENT  ra/  tlM»  R»p»H) 


Approved  for  public  release;  distribution  unlimited 


17      OlSTMiauTlON  STATEMENT  (ol  ih»  mhatfmet  aniararf  In  Block  30.   II  dlllmtt^t  Inmt  Aaport; 


H.     SU^^LEMENTAHY  NOTES 


19.     KEY  WOaOt  fConilnu*  an  fvof  »id»  II  nme*»»ir  antf  H^ntltr  ftr  ftioeft  nuMftarj 


Laminar  Flow 
Two-Phase  Flow 
Boundary  Layer 


Wedge  Flow 
Drag  Reduction 
Blowing 


20.     AaSTNACT  CCanKmM  an  rtvmtm  »ldm  It  noemmamr  mt4  HmmUtr  *f  ftlaeil  mmtbmr) 

This  thesis  theoretically  presents  the  phenomena  involved 
in  the  flow  of  an  incompressible  fluid  over  a  wedge  with  a 
second  incompressible,  lighter,  and  less  viscous  fluid  blown 
through  the  surface  of  the  wedge.   A  method  is  developed  to 
determine  the  inner  fluid  layer  thickness,  the  wall  shear 
stress  and  the  resulting  local  drag  reduction.   The  results 
predict  substantial  drag  reduction. 


OD ,  :°:*;,  1473 


eoiTlOM  O^    I  MOV  ••  IS  OaSOLETE 

S/M    0103-014- ««01    I 


UNCLASSIFIED 


ItCUHlTY  CLASai^lCATlOW  Of  Tn|«  ^AOC  (9hmn  Dmtm  Knt*T»d) 


Approved  for  public  release;  distribution  unlimited 


Two-Phase  Laminar  Boundary  Layer  Flow 
Around  a  Wedge 


by 


Dirk  van  Dord 
Lieutenant  Commander,  Royal  Netherlands  Navy 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ENGINEERING  ACOUSTICS 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  1982 


ABSTRACT 

This  thesis  theoretically  presents  the  phenomena  in- 
volved in  the  flow  of  an  incompressible  fluid  over  a  wedge 
with  a  second  incompressible,  lighter,  and  less  viscous 
fluid  blown  through  the  surface  of  the  wedge.   A  method  is 
developed  to  determine  the  inner  fluid  layer  thickness,  the 
wall  shear  stress  and  the  resulting  local  drag  reduction. 
The  results  predict  substantial  drag  reduction. 
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I.   INTRODUCTION 

The  ever  increasing  operating  depths  and  speeds  of  modern 
submarines  and  the  much  higher  resistance  to  the  effects  of 
contact  or  nearby  underwater  explosions  all  require  a  dra- 
matic increase  in  the  speed  of  underwater  weapons.   The 
existing  power  plant  systems  seem  to  have  reached  their 
limits  in  power  density  and  radical  changes  in  body  shapes 
are  limited  by  launching  facilities.   The  answer  might  lie 
in  the  field  of  boundary  layer  control.   Possibilities  are 
injection  of  polymers  into  the  boundary  layer  to  maintain 
laminar  flow,  suction  to  remove  potential  turbulence,  more 
sophisticated  weapon  shapes  to  reduce  form  drag,  or  intro- 
ducing a  gas  layer  between  the  surface  of  the  body  and  the 
liquid. 

Of  all  the  options,  the  two-phase  boundary  layer  seems 
to  be  the  most  promising  possibility  for  obtaining  an  order 
of  magnitude  increase  in  speed. 

In  this  thesis,  a  method  is  developed  to  calculate  the 
properties  of  the  laminar,  two-phase  flow  produced  by  in- 
jecting a  light,  incompressible  low-viscosity  fluid  into  the 
boundary  layer  of  the  incompressible  flow  over  a  wedge.   This 
is  just  a  step  on  the  long  way  to  a  high  speed  under^vater 
weapon. 
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II.   BACKGROUND 

A.   DISCUSSION  OF  RELATED  PAPERS 

An  extensive  literature  search  was  carried  out  and  the 
relevant  papers  are  discussed  below.   The  reduction  of  drag 
by  use  of  polymers  is  deleted  from  this  survey  because  the 
drag  reduction  expected  by  this  method  is  probably  not  suf- 
ficient to  increase  the  maximum  speed  of  the  underwater 
vehicle  substantially.   Amazingly  enough,  the  most  recent 
paper  in  the  field  of  two-fluid  laminar  boundary  layer  flow 
was  published  back  in  1969. 

1.   J.  Pretsch  [1] 

Pretsch  starts  with  Prandtl's  boundary  layer  equa- 
tions for  single-phase  flow  with  suction  or  blowing: 


dn   ^      3u    ,,3U  ^   3  u 

Ut—  +  V-r—  =  U-;—  +  V  tt 

3x     3y     jx     ,  2 

3y 


and  boundary  conditions  y  =  0:u  =  0,v  =  v;  yvRe  ^  »  :  u  =  0 
where  u  is  the  velocity  tangential  to  the  surface,  v  is  the 
velocity  normal  to  the  surface,  x  is  the  arc  length  along 
the  surface,  y  is  the  normal  to  the  surface,  U(x)  is  the  po- 
tential velocity,  v  is  kinematic  viscositv,  v  (x)  is  the 

normal  velocity  or  blowing  velocity  at  the  surface  and  Re  is 

Ux 
the  Reynolds  number.  Re  =  — . 
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For  strong  suction  (v   -^  -<»)  the  above  equations 
simplify  to  a  second  order,  ordinary  differential  equation 
of  the  form 


j2  ,                   y  V 

— sr  +  -^rr  =  0     where  r    = 

j;-2  dc,              ^ 

d^  V 


and  boundary  conditions   C  =  0  :  u  =  0 

^  -».  00  :  u  =  U  (x) 
The  solution  of  this  equation  is  the  velocity  profile 


=  1  -  e  ^ 


U(x) 


For  strong  blowing  (v  -^  °°)  the  assumption  is  made 
that  the  viscosity  term  can  be  neglected  and  the  above  equa- 
tions again  simplify  to  a  second  order,  ordinary  differential 
equation  but  now  nonlinear. 

For  a  number  of  values  for  similarity  parameters , 
Pretsch  gives  analytical  expressions  for  solutions  of  flat 
plate  and  wedge  flow. 

The  assumption  of  negligible  viscosity  is  too  re- 
strictive for  the  case  under  investigation  in  this  thesis. 
2.   E.  J.  Watson  [2] 

This  extensive  report  consists  of  three  parts,  i.e. 
Theory  of  Similar  Velocity  Distributions,  Flow  with  Uniform 
Suction,  and  Flow  with  Variation  of  Suction  Velocity.   The 
report  is  based  on  the  asymptotic  theory  of  boundary  layer 
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flow.   Although  the  author  is  aware  of  the  existence  of 
Pretsch's  work,  no  comparison  of  results  is  made. 

In  Part  I  the  author  arrives  at  a  number  of  tables 
giving  velocity  distribution  and  coefficients  of  the  series 
expansions  as  functions  of  similarity  related  constants  and 
dimensionless  parameters. 

In  Part  II  the  similarity  principle  is  abandoned  but 
the  suction  velocity  is  uniform.   The  velocity  is  expanded 
into  a  series  from  which  series  for  momentum  and  displace- 
ment thickness  and  skin  friction  are  derived.   In  the  sub- 
sequent section  a  number  of  applications  are  investigated 
and  the  results  are  totalized  in  a  number  of  tables. 

In  Part  III  the  asymptotic  theory  is  extended  to 
cover  general  two  dimensional  flow  with  arbitrary  distribu- 
tions of  mainstream  velocity  and  suction  velocity.   The 
results  are  again  displayed  in  tables.   The  author  also  de- 
rives the  results  in  Part  I  as  a  special  case  of  Part  III. 

As  this  report  only  deals  with  suction  it  is  not  of 
direct  interest  to  this  study. 

3.   R.  D.  Cess  and  E.  M.  Sparrow  [3] 

This  paper  deals  with  a  two-phase  flow  in  forced 
convection  over  a  flat  plate.   The  inner  fluid  is  obtained 
by  boiling  the  outer  fluid  at  the  wall  so  that  a  vapor  film 
is  formed.   The  authors  arrive  at  the  following  conservation 
equations: 
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f  •  '  •  +  f  f  '  '  =0 

g     g  g 

T  -  T^ 
9  '  '  +  PR  f  8  '  =  0      where  8  =  = ^ 

f  '  •  '  +  f  f  •  '  =0 

Li  Li   Li 

with  boundary  conditions: 

plate  surface:  f  (0)=f  '(0)  =0,  9(0)  =1 

g       g 


interface:     f^(0)  =  ^f^  j  fglng^) 

Li 

(py)  '"^ 


e(nj)  =  0 

free  stream:  f  '  '  ^  2  as  n^  -^  <^ 

where  f   and  f^  are  the  stream  functions  for  respectively 
g      L 

inner  and  outer  fluid,  T  is  the  local  temperature,  T,  is  the 

w 

temperature  at  the  wall,  T^  is  the  temperature  in  the  free 

stream,  PR  is  the  Prandtl  number,  o  is  the  mass  density  and 

U  is  the  dynamic  viscosity.   f   is  a  function  of  the  dimen- 

sionless  parameter  n   defined  by  n   =  ^  / .   n  •  is  the 

^  g  ^   g   2^JVgX     'gi 

value  for  this  parameter  at  the  interface  of  the  two  fluids. 
f_  is  a  function  of  the  dimensionless  parameter  n,  defined 

/IT"  ^ 

by  n.,  ~  2  /  — X  *   "^^^  subscripts  L  and  g  refer  to  respective- 
ly outer  and  inner  fluid.   The  primes  denote  differentiation 


with  resDect  to  the  relevant  dimensionless  parameter.   This 
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problem  is  solved  by  using  an  Emmons-Leigh  [4]  table  and  a 
successive  approximation  procedure. 

4.  K.  Nickel  [5] 

In  this  paper,  the  laminar  flow  of  an  incompressible 
fluid  without  heat  transfer  is  investigated  using  a  method 
of  bounding  the  solutions  applied  to  Prandtl's  equations. 
After  introducing  the  well-known  dimensionless  parameters  n 
and  f(n)/  Prandtl's  partial  differential  equations  reduce  to 
the  ordinary,  third  order  non-linear  differential  equation: 

f*"  +  ff'  +  3(1  -  f'^)  =  0 

with  boundary  conditions:   f (0)  =  f^ 

f '  (0)  =  0 

f  (oo)  =  1 
The  solutions  are  reflected  in  a  number  of  graphs, 
especially  "abbildung  8"  giving  f''(0)  as  a  function  of  3 
and  f|-.  will  prove  to  be  very  useful  in  this  thesis. 

5.  R.  D.  Cess  and  S.  M.  Sparrow  [6] 

This  paper  is  an  extension  of  the  paper  discussed  in 
Section  3  above  and  deals  with  the  case  of  a  subcooled  fluid. 
To  extend  the  earlier  results  it  was  necessary  to  consider 
the  conservation  of  energy  within  the  fluid.   The  results  are 
again  reflected  in  a  number  of  figures . 

6.  E.  M.  Sparrow,  V.  K.  Jonsson  and  S.  R.  G.  Eckert  [7] 
In  this  paper  the  laminar  tv/o-phase  flow,  obtained  by 

injecting  gas  into  a  liquid  at  the  wall  of  a  flat  plate,  was 
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investigated.   The  authors  arrive  at  the  following  ordinary 
differential  equations 
for  the  gas:        ^a'"  "^  ^a^a' '    =  ^ 
for  the  liquid:     ^ l  '  '  '  "*"  ^L^L*'  ^  ^ 
with  boundary  conditions: 

at  the  wall:        ^cr^^^  ~  ^aO    ^^^^^^^^   parameter) 

f  ' (0)  =0   (no  slip  condition) 

at  the  interface:   f  (n  •)  =  f-r  (0)  =  0  (no  mass  transfer) 

g   gi     L 

f  ''(n  .)  =Rf^''(0)  (stress  continuous) 
g    gi       L 


where  R  = 


f  ' (n  ■)  =  f X ' (0)   (velocity  continuous) 
g    gi      L 

for  n,  ^  =°:        f  '  (n-r)  "^  2   (velocity  approaches  U  ) 

The  authors  solve  these  equations  by  applying  a 
series  expansion  and  subsequently  a  successive  approximation 
procedure,  similar  to  the  one  described  in  the  paper  by  Cess 
and  Sparrow  (Section  3  above) . 

The  results  are  given  in  a  table  and  a  graph.   It  is 
proven  that  evaporation  at  the  interface  only  works 
advantageous . 

7.   H.  P.  Eichenberger  and  J.  D.  Qffutt  [3] 

Under  assumption  of 

-  shear  stress  in  liquid  is  negligible 

-  interfacial  velocity  and  pressure  distribution  are  those 
imposed  by  the  potential  flow 
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-  the  inertia  of  the  gas  is  negligible 

-  the  gas  flow  is  laminar 

-  body  moves  steadily  through  water  at  zero  incidence 

-  gas  and  water  properties  are  constant;  in  particular  the 
density  of  the  gas  remains  constant 

the  problem  of  a  two-phase  flow  over  different  kinds  of 

bodies  is  solved  by  a  finite  differences  method.   Results 

are  given  in  a  number  of  graphs . 

8.   W.  S.  Bradfield,  R.  0.  Barkdoll  and  J.  T.  Byrne  [9] 

This  paper  discusses  experiments  with  a  "gassing" 

model  to  investigate  possible  drag  reductions.   The  vapor 

layers  were  obtained  by  film  boiling,  sublimation  of  a  solid 

surface  and  chemical  reactions.   The  experiments  indicated 

that  for  film  boiling  a  region  exists  for  which  a  stable 

vapor  layer  may  be  expected. 

Assuming 

-  steady,  two  dimensional  laminar  flow 

-  pressure  everywhere  constant 

-  heated  surface  temperature  constant 

-  temperature  of  liquid  vapor  interface  is  equal  to  satu- 
rated liquid  temperature 

-  thermodynamic  and  transport  properties  of  liquid  and 
vapor  are  constant 

-  thermal  radiation  effect  on  vapor  and  liquid  layers  is 
neglected 

-  buoyancy  effect  on  forced  convection  boundary  layers  is 
neglected 

-  velocity  and  temperature  profiles  are  linear  in  vapor  and 
parabolic  in  liquid 
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and  using  the  Pohlhausen  integral  method  the  results  of  the 
experiments  are  predicted.  In  some  of  the  experiments  dra- 
matic instabilities  were  found,  and  the  authors  had  to  con- 
clude that  producing  a  stable  vapor  layer  between  a  surface 
and  a  liquid  drastically  reduces  friction  drag  but  unstable, 
two-phase  flow  (like  nucleate  boiling  and  uncontrolled  sub- 
limation) may  actually  cause  the  drag  to  increase. 
9.   T.  J.  Lang  [10] 

In  this  patent  Lang  "designed"  a  torpedo  capable  of 
using  the  gas  film  drag  reduction  phenomena.   The  lift  is 
provided  by  the  body  shape  and  the  exhaust  gases  are  used  to 
produce  the  gas  layer.   At  the  rear  end  of  the  torpedo  the 
gas  layer  is  removed  by  suction. 
10.   E.  R.  van  Driest  [11] 

The  author  indicates  that  the  quest  for  higher 
speeds  of  underwater  craft  can  be  solved  by  using  the  knowl- 
edge acquired  in  high  speed  aerodynamics. 

In  a  number  of  figures  the  author  gives  local  fric- 
tion coefficients  as  a  function  of  Reynolds  number,  speed, 
wall  temperature  and  ocean  depth.   Further  effects  of  hydro- 
dynamic  heating,  cavitation  and  insertion  of  additives  are 
discussed. 

The  author  concludes  that  although  the  propulsion 
systems  have  to  be  developed,  it  is,  in  principle,  possible 
to  obtain  high  velocities  for  underwater  vehicles. 
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B.   QUO  VADIS 

Considering  the  above  discussion  it  was  felt  that  more 
knowledge  had  to  be  obtained  of  the  phenomena  taking  place 
in  the  two-phase  boundary  layer.   A  first  step  would  be  to 
extend  the  work  on  flat  plates  to  two-dimensional  bodies 
with  pressure-gradients.   It  is  logical  to  follow  the  paper 
of  Sparrow,  Jonsson  and  Eckert  [7] ,  discussed  above  in  Sec- 
tion A-6. 
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III.   THEORY 


A.   GOVERNING  EQUATIONS 

To  predict  the  skin  friction  and  drag  characterics  from 
first  principles  the  starting  position  is  based  on  the  basic 
conservation  laws  governing  the  velocity  distribution,  here 
in  Cartesian  coordinates/ 


Continuity; 
Navier-Stokes : 


at     dX       3y 


Du  _      5P    3 

Dt        3x    3x 


aZ 


I  t3u   2  , . 
y  2^-  -  ^div  w 
[  3x   3 


i] 


^  3  r  f5u  ^  3vll  ^  3  r  f3w  ^  3ul 
Dt        ^y    c-y|^^[  3y    3      J  J 


3z 


, 3v  ^  3w 


^  3    f3u  ^  3vl 


D^  -  -7    3?  ^  3 
^Dt  -  ^  -  T^  ^  TI 


y  .  2t —  -  ^div  w 
jz    3 


3X 


3w  ,  3ul 


+ 


I  JV    J^ 

^I3z    3yj 


Eqn  of  State; 


P  =  p  R   T 
U  =  U(P,T) 


do    ^^T    dW 
Energy  Equation;  -^   =  -rrr^  + 


dt 


dt 


dt 


where  w  is  the  velocity  with  components  u,  v,  w  in  respec- 
tively the  X,  y  and  z  direction,  —  is  the  operator 

!^  '^  Ci  '^ 

■j£   +  u-^  +  v^  +  w-r^    ,    X,    Y   and  Z  are  the  body  forces  in 
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respectively  x,  y  and  z  directions,  P  is  the  local  pressure, 
R  is  the  gas  constant,  Q  is  the  heat  into  the  system,  E-,  is 
the  internal  energy  and  W  is  the  work  done  by  the  system. 

Assuming  steady  state,  absence  of  body  forces,  viscosity 

^P 
u(T,P)  is  constant,  two-dimensional  flow,  -—-   is  negligible, 

ay 

laminar  flow,  and  both  fluids  incompressible  (the  energy 
equation  decouples) ,  these  equations  simplify  to  the  Prandtl 
boundary  layer  equations: 

Continuity:     -|H  +  -^  =  0 

dx    dy 

2 

Navier-Stokes:   u^H  +  v-^  =  -  ^  -p  +  v^ 

oX     iy      0  3x     -2 
^  3y 

Using  the  subscript  g  for  properties  of  the  inner  fluid  and 

the  subscript  L  for  the  outer  fluid,  it  is  possible  to  write 

two  sets  of  equations  describing  the  two-fluid  laminar 

boundary  layer. 

Inner  fluid:   iH  +  4^  =  0  (1) 

c*x    3y 

u|»  +  v4H  =  -  i  4£  +  V  £^  (2) 

9x     ay      p  ax     g^  2  ' 

Outer  fluid:   iH  +  ^  =  0  (3) 

dx         5y  ^    ' 

5u  _,  au    1  ap  ^   3^u 

Urr—    +    V-r-  =  -  —  ^  +  v_ s-  4) 

3x     oy      p  cx     L^  2 
^  3y 

To  describe  the  problem  fully,  however,  it  is  required  to 

state  the  boundary  conditions  and  the  conditions  at  the 

interface  between  the  two  fluids.   At  the  wall,  in  the  inner 
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fluid,  the  velocity  parallel  to  the  wall,  u,  is  zero  (no 
slip)  and  the  velocity  perpendicular  to  the  wall  is  equal  to 
the  imposed  blowing  velocity.   At  the  interface  the  tangen- 
tial velocity  and  shear  stress  have  to  be  continuous.   Fur- 
thermore, the  mass  transfer  across  the  interface  is  assumed 
to  be  zero.   In  the  outer  fluid  u  has  to  approach  U(x)  as  y 
approaches  infinity.   A  formal  statement  of  these  conditions 
is: 
at  y  =  0        u  =  0 


V  =  V 


w 


at  interface    u   =  u_ 

g    L 


3u. 


fSu, 


gl  3y  Ji   Li  3y  J^ 

m  =  0  or  V   (interface)  =  v   (interface)  =  0 

g  L 

in  outer  fluid   u  -*  U(x)  as  y  ->  =0 

For  the  given  assumptions,  equations  1-4  and  the  boundary 
conditions  describe  the  problem  completely.   This  system  pre- 
sents substantial  difficulties  in  attempting  a  general  solu- 
tion.  Similarity  solutions,  however,  will  solve  the  deadlock 

B.   TRANSFORMED  EQUATIONS 

Similarity-type  solutions  are  well  known  and  have  the 
characteristic  that  the  velocity  profiles  have  similar  shapes 
at  all  x-positions  along  the  surface.   Mathematically  this 
gives  the  enormous  simplification  that  the  governing  partial 
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differential  equations  are  converted  into  ordinary  differen- 
tial equations.   Before  defining  the  dimension less  para- 
meters it  should  be  noted  that  the  continuity  equations  are 
satisfied  by  the  stream-function  '4;  defined  by 


u  =  -r-^  and  V  =  -  -r^ 
3x  jX 


Following  Schlichting  [12]  ,  the  following  dimensionless 
parameters  can  be  introduced  for  the  inner  fluid: 

Reynolds  number   Re   =  (where  U   and  L  are  some 

^        g 


reference  quantities) 


n   =  ■= — r—t        t,      =  zr      (where  g(x)  is  a  dimensionless  scale 
'g   Lg  (x)    ^g   L         -a  ^  ' 

factor  for  y) 


'J;(x,y)  yRe 
f  (n  ,C  )  =  ^  ,,,  ,   ,  ?  or  iJ(x,Y)    =   f  (n  ,>  )L  U(x)  g(x) 
'g'^g'    L  U(x)  g(x)       ^'  '^'  g   g     /^o^- 

/   g 

To  simplify  further  manipulation,  it  is  wise  to  state  at 

this  point  that  to  ensure  similarity  f(n  ,  E,    )    can  not  depend 

g  g 

on  5  =   E,    (x)  so  that  f  =  f  ( n  )  . 
^g    ^g  g 

From  the  above,  the  velocity  components  may  be  derived, 
u  =  U(x)  f '  (rig) 

and  V  =   -Li^^(U(x)g(x))    ^  hJll2SL^l2Sl  ,    f    .(n    )Ag(^) 

/   g  *^   g 

Inserting  these  expressions  into  the  Prandtl  boundary  layer 
equation  and  after  a  considerable  amount  of  algebra  the  equa- 
tion converts  into 
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where  a  =  ^-^i^   ^^^^^^  '^^^^^ 

00 

o   -   L  g^(x)  dU(x) 

'^   ~        U      dx 

00 

Reasoning  in  the  same  way  for  the  outer  fluid  after 
defining 

U  L 


=  (y  -  s) 


X 


^®L  ~  ~7r    '  "^L  "  L  g(x)  /^^L  '   ^L    L 


fr  (Hr  ,Ct)   = 


•L"'L'*"L'  ~  L  U(x)  g(x) 
one  arrives  at 


where  a  and  3  are  defined  the  same  as  for  the  inner  fluid. 

It  should  be  noted  that  to  satisfy  the  similarity  re- 
quirement  a  and  3  should  not  depend  on  x.   Still  following 
Schlichting  [12],  it  is  possible  to  arrive  at  the  following 
general  conditions  for  U(x)  and  g(x): 

8 


2   - 


U(x)  ^  K2a-3|(2a  -  3)? 


r-^   ^'L 


2  a- 8 

where  K  is  an  integration  constant 


g(x)  =  _y(2a  -  3)f  ^^ 


Now,  without  loss  of  generality,  take  1=1  and  introduce  a 


new  constant  m  =  ■=— ^  so  that  the  expressions  for  —r, — -   and 


UU) 

u 


00 


g(x)  convert  to 
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U(x)  ^  ^(1+m)  |'_2_  xl'^ 
U  1+m  lJ 


2   X 


"^^^^  "  \jm+l    L  U(x) 
Furthermore,  the  transformation  equations  for  the  ordinates 


/m+1  U(x)     ^       ,      j:^     /m+1  U(x) 
are  n  =  y  ~  /-  -^  / 


and  n-j^  =  (y-o)  /— 2" 


•g    ^y  2   VgX        L     -    V  ^   ''l^ 
With  Schlichting  [12],  it  can  be  noted  that  these  re- 
quirements can  be  satisfied  by  the  flow  in  the  neighbourhood 
of  the  stagnation  point  of  a  wedge  with  included  angle  3tt, 
where  U(x)  is  given  by  Cx  . 

Next  the  boundary  and  interface  conditions  have  to  be 
converted  by  expressing  them  into  the  newly  defined  dimen- 
sionless  parameters 
At  y  =  0 

or  n  =  0     u  =  U(x)  f  '  (n)  =  0   or  f  '  (0)  =  0 
y  y  y 

-L  f  (n  )  , 

vRe 


-v   'Re_ 

r   /    \  W  y'     g 

or  f_{n„)  =  ^- 


g   g     L  d(U(x)  g(x)  ) 

dx 

at  interface  n   =n-   m  =  Oorf(ri-)=f-r(0)=0 
g    gi  g   gi     L 


n^  =  0 


"g  =  ^L  °''   fg''^gi>  =  ^l''°' 


du       du 

^g  ^  =  ^L  ^7  ^^  ^g'  '  <^i)  =  ^  "'l'  '  ^°^  ^^^^^^  ^  = 

for  y-^^^orri^-oo  u->U  (x)  or  f  t.  '  ->  1  as  r\-  -^  ^ . 

Li  Ij  L 
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It  should  be  noted  that 

-  to  ensure  that  f  (n  )  is  independent  of  x: 

y  y 

y     o.  d(U(x)  g(x)  )  ^    Inserting  the  expressions  for  U(x) 
w        dx  ^        - 


and  g(x)  one  obtains  v  "^^  J- 

w   V 


U{x) 


-  R  depends  on  the  properties  of  both  fluids  depending  on 
temperature  and  pressure. 

Summarizing  the  problem  can  be  formally  stated  by 

equations:   f'''+f   f''+3(l-f'^)=0  (5a) 

g        g   g  g 

and  conditions:   f  ' (0)  =0 

g 


f^(0)  =  -V,  !  ,„,,,^^  ,,,    ,  where  v  '^>  '^ ''^^ 


■g'  w  '  (m+1)  V  U(x)  ™^^=^^  ^w    V  x 
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IV.   EXPLORED  WAYS  OF  SOLUTION 

A.   SERIES  SOLUTION 

Following  the  procedure  of  Sparrow,  Jonsson  and  Eckert 
[7],  a  series  solution  of  the  momentum  equation  seems  to  be 
appropriate.   Expanding  f  (n  )  with  respect  to  n   =0  gives 

f  '  '  (0)n  ^   f  •  •  '  (0)n  ^ 

The  first  two  terms  in  the  expansion  are  available  from  the 
boundary  conditions.   The  f  '' (0)  is  unknown  for  now  and  will 
be  determined  by  the  solution.   Higher  derivatives  can  be 
expressed  in  terms  of  lower  derivatives  by  taking  sequen- 
tially derivatives  of  both  sides  of  the  momentum  equation 

f • • '  =  -f  f • •  -  3(1  _  f .2) 


and  evaluating  at  i   =  0 . 

Some  expressions  for  derivatives  are: 

f <^^^)  (0)  =  -f (0)f ' '  (0)  -  3 

f (^^^  (0)  =  f^(0)f ' ' (0)  +  3f (0) 

f^^No)  =  -f^(0)f"(0)  -  3f^(0)  +  (23-l)f '  '^(0) 

f(^^^{0)  =  f'^(0)f"(0)  +  3f^(0)  -  (3S-5)f  (0)f  •  '^(0) 

-  (63-4) 3f '  '  (0) 
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It  was  impossible  to  find  a  recurrency  in  the  expressions. 
On  applying  the  interface  condition  f (n  ■)  =  0  makes  it  pos- 
sible  to  write 


f  "(0)  =       -^^^^ 


'^gi    L^      f  '  '  (0)k! 
k=3    ^ 


It  is  now  possible  to  express  f  ' ' (0)  in  liquid  and  gas  para- 

y 

meters.   Here  Sparrow,  Jonsson  and  Eckert  refer  to  an  earlier 
paper  by  Cess  and  Sparrow  [3] .   This  paper  is  discussed  in 
Chapter  II,  Section  A3.  This  paper  gives  a  table  relating 
f_  '  (0)  and  f^ ' '  (0)  which  corresponds  to  a  solution  of  the 

Li  Li 

equation  for  a  flat  plat. 

At  this  point  it  was  realized  that  the  data  necessary  to 
go  any  further  was  not  available,  therefore  this  method  was 
abandoned,  and  other  methods  of  solution  were  explored. 

B.   NUMERICAL  SOLUTIONS 
'  1.   lODE  [13] 

lODE  is  an  interactive  ordinary  differential  equation 
solver  developed  by  the  Computer  Center  of  the  Naval  Post- 
graduate School.   The  program  accepts  a  number  of  first  order 
non-linear  ordinary  differential  equations  with  associated 
initial  conditions  and  solves  the  equations. 

In  this  case  equations  5  were  converted  into: 

dy^ 

-^  =  y2  =  f    Note:  Y^^   ^  (6a) 
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dy 


2 


dn    ^3 
dy. 


=  y,  =  f "  (6b) 


2 


__  =  f...   =  -y^y^  -   3  +  6y2  (6C) 

where  y,  ,  Yj    ^Jid  y^  are  duittmy  variables  for  the  derivatives 
of  f{r\),    with  initial  conditions 

Y^i^)    =  f(0)  (7a) 

y2(0)  =  f  (0)  (7b) 

y3(0)  =  f  '  (0)  (7c) 

The  equations  6  are  the  same  for  both  fluids  but  as  noted 

before  the  initial  conditions  differ. 

Sanders  [14]  suggested  that  it  should  be  possible  to 

guess  the  thickness  of  the  inner  fluid  layer  n  .  and  then  a 

value  of  f  ' ' (0)  could  be  found  by  trial  and  error  so  that 

no  mass  transfer  occurs  across  the  interface  (f  (r,  .)  =0). 

g   gi 

The  values  for  f  ,  f   and  f'   at  the  interface  give  the 

g   g      g  ^ 

initial  conditions  for  the  outer  fluid  equations  6 .   The 
values  of  f_'  (n^)  should  approach  unity  for  large  enough  r\ 

ij      Ij  Li 

while  f T  '  '  (n-p)  approaches  0.   If  f  '  (n,.)  does  not  go  to  unity 

Li        Li  Li      Li 

the  initial  guess  for  the  layer  thickness  was  wrong  and 

another  value  for  the  inner  fluid  layer  thickness  -n    is 

^  gi 

chosen.   This  process  is  repeated  until  a  layer  thickness  is 
found  for  which  f^  '  (°°)  =1  and  f_  '  *  (=")  =0. 

Li  Li 

The  calculations  carried  out  on  the  IBM  3033A?  com- 
puter system  showed  that  this  procedure  worked,  but  while  the 
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calculations  were  fast  the  interactive  side  of  the  program 
has  proven  to  be  too  time  consuming. 
2 .   Desk  Top  Computer 

So  the  necessity  arose  to  find  some  way  to  decrease 
the  time  spent  communicating  with  the  computer  and  to  auto- 
mate the  trial  and  error  process.   Furthermore,  it  was 
necessary  to  find  some  way  to  interrupt  the  execution  and 
to  execute  the  program  by  steps  to  determine  instabilities 
in  the  solution. 

Since  a  Wang  2200  desk-top  computer  with  line  printer 
was  available,  the  problem  was  rewritten  for  this  computer 
and  the  solution  for  a  number  of  cases  was  carried  out.   It 
has  to  be  realized  that  a  desk-top  computer  is  relatively 
slow  so  that  a  compromise  between  calculation  time  and  ac- 
curacy of  the  result  had  to  be  found.   It  was  decided  that 
the  primary  goal  was  to  find  a  method  of  solution  and  not  to 
tabulate  a  large  amount  of  accurate  data.   In  the  light  of 
all  the  assumptions  high  accuracy  would  be  doubtful  anyway. 
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V.   SOLUTIONS 

A.   DEFINITION  OF  THE  PROBLEM 

The  momentum  equation  will  be  written  as  three  first 
order  differential  equations  as  for  the  "lODE  program": 

dy. 


^  =  y.  =  f 


dn    -^2 

dy^ 

=  y,  =  f" 


dn 
dy 


3  0    o   2 

=  -yiy^  -  3  -  3y. 


dn     n^3    "   "^1 

Now  multiplying  both  sides  of  these  equations  by  dn  and 
changing  to  finite  differences  one  arrives  at 


Ay^  =  y-,  -An 


2    ^3 

2 
Ay3  =  (-yi*y3  -  3  -  3y^  ) -An 

Now  by  stepping  n  through  the  desired  range  of  values, 
starting  with  known  or  guessed  initial  values  for  '^  /  y-,  f  y2 
and  y-.  it  is  possible  to  calculate  y,  ,  y-  and  y^  as  a  func- 
tion of  n- 

It  has  to  be  noted  at  this  point  that  we  approximate  the 
integrals  by  a  "Lower  Riemann  Sum".   It  is  known,  however, 
I   that  using  a  "sufficient"  number  of  steps  guarantees  a  satis- 
factory accuracy.   This  procedure  is  applicable  for  both  the 
inner  and  the  outer  fluid.   The  initial  values  for  the  outer 
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fluid  follow  logically  from  the  solution  of  the  inner  fluid 
and  the  interface  conditions. 

Given  values  for  the  blowing  parameter  f  (0)  and  the 
properties  of  the  fluids  parameter  R  the  procedure  is  as 
follows: 

-  Choose  an  inner  fluid  layer  thickness. 

-  Find  f'  (0)  by  an  interval  halving  method  so  that  f  (r,  .) 

(no  mass  transfer  across  the  interface)  equals  0 ,  or  "^ 
in  our  program  is  smaller  than  a  chosen  value. 

Note:   f(0)  =  given;  f'(0)  =  0  and  f''(0)  is  chosen  as 
indicated. 

-  Apply  the  interface  conditions  on  the  values  for  t,    f 
and  f'  at  the  interface  to  find  the  initial  values  for 
the  solution  of  the  outer  fluid  equations.   Choose  an 
end  value  for  r\      large  enough  to  be  well  "outside"  the 
boundary  layer. 

-  Check  if  f^'  (^)  =  1  but  also  if  f  '  '  (=°)  =  0  as  it  is 
possible  that  an  unstable  solution  is  found.   In  case 
the  conditions  are  satisfied  the  solution  is  found.   If 
not  another  guess  on  the  inner  fluid  layer  thickness  has 
to  be  made. 

-  The  value  of  the  inner  fluid  layer  thickness  is  chosen 
in  another  interval  halving  process. 

-  This  process  is  repeated  until  the  requirements  of 
f_  '  (^)  =  1  and  f_  '  '  (°°)  =0  are  met. 

Li  Li 

A  copy  of  this  program  is  enclosed  as  Appendix  A. 

From  the  above,  it  follows  that  the  initial  intervals 

for  f  ' ' (0)  and  n  •  have  to  be  determined.   The  interval  for 
g  gi 

n  .  should  be  as  narrow  as  possible  and  the  interval  for 
f'  (0)  should  be  a  little  wider  than  expected  to  prevent  in- 
stabilities and  "limit-riding".   These  values  were  determined 
by  trial  and  error  and  the  results  are  discussed  in  a  later 
chapter. 
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B.   TESTING  THE  PROGRAM 

1 .  Number  of  Steps  of  Independent  Variable 

The  program  is  used  with  100  steps  of  the  independent 
variable.   A  sample  problem  was  solved  with  both  100  and  200 
steps.   The  deviations  were  in  the  fourth  decimal  (see  Table 
5.1).   The  time  required  for  the  calculation  however  was 
doubled  from  approximately  45  minutes  to  90  minutes. 

2.  One  Fluid  Flow 

It  is  clear  that  when  R  approaches  unity,  the  two 
fluids  are  hydrodynamically  identical.   Computer  solutions 
were  obtained  for  a  single  fluid  flow  for  3=0  (flat  plate) 
and  B  =  1  (stagnation  flow)  and  further  for  3  =  0.2,  0.4, 
0.6  and  0.8.   The  results  were  in  good  agreement  with  Nickel 
[5],  "abbildung  8".   (See  Table  5.2.) 

3 .  Solution  for  R  Close  to  1 

A  solution  with  a  R  of  1.1  was  compared  to  that  for 
a  one-fluid  problem  and  the  results  were  close  enough  to 
give  further  confidence  in  the  correctness  of  the  program  to 
handle  non-unity  values  of  R.   (See  Table  5.3.) 

4 .  Solution  for  Small  Wedge  Angles 

The  solution  for  a  3  of  0.01  was  close  enough  to  that 
of  the  flat  plate  to  give  confidence  about  the  capabilities 
of  the  program  to  correctly  handle  non-zero  values  of  S. 
(See  Table  5.4.) 
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VI.   RESULTS 

A.  DEFINITION  OF  A  SAMPLE  PROBLEM 

To  demonstrate  the  capabilities  of  the  program,  a  special 
case  was  chosen  for  thorough  investigation.   For  this  case, 
"air  at  standard  conditions"  and  water  were  chosen  as  the 
two  fluids  so  that  R  =  200.   For  a  range  of  blowing  para- 
meters f  (0),  and  wedge  angles  3,  the  velocity  profile  was 
y 

calculated  and  the  layer  thickness  n  .  of  the  inner  fluid 

go. 

and  the  shear  stress  parameter  f  ' ' (0)  at  the  wall  were 

g 

determined. 

B.  THICKNESS  OF  THE  INNER  FLUID  LAYER 

Figure  6 . 1  gives  the  dimensionless  thickness  of  the  inner 
fluid  layer  n  •  as  a  function  of  the  blowing  parameter  f  (0) 
for  various  wedge  angles  3.   Figure  6.2  is  an  enlargement  of 
these  curves  for  small  blowing-parameters.   These  graphs  lead 
to  the  following  observations : 

-  For  small  values  of  the  blowing  parameter  the  curves  are 
straight  lines,  ngj_  decreases  very  little  with  increasing 
wedge  angle  and  is  directly  proportional  to  f  (0). 

-  For  larger  blowing  parameters  fg(0),  nqj_  increases  more 
rapidly  with  increasing  fg(0).   For  large  blowing  para- 
meters and  small  angles  tne  thickness  of  the  inner  fluid 
layer  increases  drastically  with  increased  blowing 
parameter. 

-  For  large  values  of  the  blowing  parameter,  the  thickness 
of  the  inner  fluid  layer  decreases  rapidly  with  increas- 
ing wedge  angle  3. 
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C.   WALL  SHEAR  STRESS 

Figures  6.3  and  6.4  give  the  dimensionless  wall  shear 
stress  f  ' ' (0)  as  a  function  of  the  blowing  parameter  f (0) 

g 

and  the  wedge  angle  3 . 

-  For  all  values  of  wedge  angle  S,  blowing  decreases  the 
wall  shear  stress  fg'MO).   This  reduction  is  spectacular 
for  small  values  of  the  blowing  parameter;  e.g.  for  both 
the  flat  plate  (3  =  0)  and  stagnation  flow  (3=1) 

fg' ' (0)  is  reduced  to  half  the  non-blowing  value  when 
fg(0)  =  0.0025. 

-  For  all  values  of  the  blowing  parameter  the  wall  shear 
stress  f'  (0)  increases  with  increasing  wedge  angle  3; 
the  smaller  thickness  of  the  inner  fluid  layer  caused  by 
the  higher  potential  velocity. 

-  For  large  values  of  the  blowing  parameter,  the  wall  shear 
stress  is  very  small  and  nearly  equal  for  all  wedge 
angles.   The  values  obtained  indicate  that  the  wall  shear 
stress  for  large  blowing  parameters  approaches  the  wall 
shear  stress  for  a  wedge  in  an  infinite  volume  of  inner 
fluid.   (See  Table  6.1.) 


D.   VELOCITY  PROFILES 

Figure  6.5  gives  velocity  profiles  (f'(n)  vs.  n)  for 

various  blowing  parameters  for  a  wedge  angle  3  =  0.2  (36'^). 

The  r\-    and  n   have  been  expressed  in  the  same  relative  scale 
L     g       ^  

y/Re  (y-o)^'Re^ 

by  noting  that  ^^  =  j;-^     ^^d  n^  =   ^  g(^)    resulting  in 

An^  =  3  • 

The  following  observations  can  be  made: 

-  For  the  smaller  blowing  parameters  the  velocity  profiles 
in  the  inner  fluid  are  nearly  straight  lines. 

-  The  stronger  the  blowing  the  thicker  the  layer  of  the 
inner  fluid  causing  a  larger  proportion  of  the  potential 
velo"city  to  be  obtained  in  the  inner  fluid  layer.   For 
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small  blowing  parameters  this  proportion  of  the  potential 
velocity  increases  rapidly  v;ith  increasing  blowing  para- 
meter.  Figure  6.6  gives  the  proportion  of  the  potential 
velocity  (or  f  ' (n  •))  as  a  function  of  the  blowing  para- 
meter f  (0)  .   "^   ^^ 

g 

Particularly  interesting  is  the  fact  that  the  velocity 
profiles  in  the  inner  fluid  and  the  outer  fluid  have  op- 
posite curvatures.   This  is  comparable  to  the  boundary 
layers  associated  with  the  gas-filled  cavity  attached  to 
a  body  after  water  entry.   The  cavity  around  such  a  body 
has  a  stable  interface;  breakdown  occurs  only  when  too 
much  gas  is  left  at  the  rear  separation.   Although  no 
rigorous  proof  can  be  given  it  is  felt  that  this  charac- 
teristic of  the  velocity  profile  provides  stability  of 
the  interface  in  two  fluid  boundary  layer  flow. 


E.   DRAG  REDUCTION 

To  characterize  the  amount  of  drag  reduction,  the  local 
drag  parameter  f  ' ' (0)  obtained  for  the  two  fluid  flow  was 
compared  to  the  local  drag  parameter  obtained  without 
blowing. 

To  compare  these  local  drags,  a  minor  converson  has  to 
be  made.   The  shear  stress  at  the  wall  is 

In  terms  of  the  transformed  variables,  the  expression  for  t 
becomes 

U(x) 


which  can  be  expressed  in  the  outer  fluid  parameters 


U(x) 


L  L  g(x)  y   L  I  (py)^  I    g 


h 
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For  the  outer  fluid,  without  blov/ing,  the  reference  shear 


stress  Tq  is: 


-0  =^lZ^/^l   ^^l"  (^)^=o 


f  (0)=0 

g 


Now  taking  the  ratio: 


^0    ^(fL"<^)^f  (0)  =  0 

g 

The  values  of  this  reference  shear  stress  parameter  f  ' ' (0) 

are  calculated  and  reflected  in  Figure  5.7.   Since  R  is  a 

known  constant,  for  each  value  of  the  shear  stress  parameter 

f  '• (0)  it  is  possible  to  determine  the  ratio  t/t^.   Table 
g  u 

6.2  gives  the  values  of  t/t^  as  a  function  of  wedge  angle  S 
and  blowing  parameter  f  (0)  for  R  =  200. 

Inspection  of  this  table  reveals  that  a  drag  reduction 
is  observed  for  all  blowing  param.eters  and  wedge  angles.   It 
is  remarkable  that  for  small  values  of  the  blowing  parameter 
the  ratio  t/t.  decreases  for  increasing  wedge  angle  3  while 
for  large  values  of  the  blowing  parameter  f (0)  the  ratio 
t/t^  increases  with  increasing  wedge  angle  S.   For  all  wedge 
angles  3/  however,  the  ratio  t/Tq  decreases  with  increasing 
blowing  strength. 


F.   REQUIRED  POWER 

The  power  required  to  maintain  blowing  is 

Power  =  P  •  V 
w 

where  P  is  the  local  pressure  given  by 
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^     STAG   [jxJ 

where  P_^^_  is  the  stagnation  pressure  at  the  leading  edge 
STAG 

of  the  wedge  and  ||  =  - p  u||  =-p  C^ m  x^^"^  <  0.   If  the 

inner  fluid  is  at  stagnation  pressure  P^_^-  a  pressure  dif- 

3P 
ference  of  ^—  at  the  wall  exists  so  that  the  fluid  flows  out 

oX 

by  itself.   The  flow  rate  then  can  be  controlled  by  the 
blowing  hole's  size.   Of  course,  the  ejected  fluid  has  to  be 
replaced  requiring  it  to  be  brought  up  to  P__^^.   In  practice 

o  i  AG 

this  is  of  less  interest  as  a  considerable  amount  of  exhaust 
gases  of  the  propulsion  systems  at  high  pressures  should  be 
available . 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  shown  theoretically  that  a  dramatic  drag 
reduction  is  possible  for  the  infinite  wedge  by  blowing  with 
a  lighter  and  less  viscous  fluid,  and  that  the  parameters 
involved  can  be  calculated  with  relatively  simple  means. 
Furthermore  the  developed  method  gives  another  way  of  solving 
the  flat  plate  case  and  the  stagnation  flow  case. 

This  work  directly  invites  a  number  of  follow-on  investi- 
gations of  phenomena  associated  with: 

-  finite  two-dimensional  bodies 

-  three-dimensional  bodies 

-  a  compressible  inner  fluid 

-  experimental  verification  of  the  theory 

-  stability  of  the  interface 
all  of  equal  importance. 
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APPENDIX  A 

"BASIC"  PROGRAM  FOR  WANG  DESK  TOP  COMPUTER 
TO  SOLVE  TWO  FLUID  LAi4INAR  FLOW  CHARACTERISTICS 


1  REM  THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  F,F',F" 

2  REM  AND  NGI  FOR  A  TWO-LIQUID-LAMINAR-BOU^JDARY-LAYER 

3  REM  GIVEN  A  BLOWING -PARAMETER  F(0),A  WEDGEANGLE  B  AND 

4  REM  A  LIQUID-PROPERTIES-CQNSTANT  R. 

5  REM 

6  REM  TO  USE  Tl-C  PROGRAM  ONE  HAS  TO  GUESS  ON  AN  INITIAL 

7  REM  INTERVAL  FOR  THE  INNER-FLUID-LAYER-THICKNESS  AND 

8  REM  FOR  THE  LOCAL  SHEARSTRESS-PARA^CTER  F"(0). THESE 

S    REM  VALUES  CAN  BE  OBTAINED  FROM  GRAPHS  IN  THE  THESIS. 

10  REM  ONE  HAS  TO  TAKE  THE  INTERVAL  FOR  THE  THICKNESS  AS 

11  REM  SMALL  AS  POSSIBLE  AND  THE  INTERVAL  FDR  THE  SHEAR - 

12  REM  STRESS  A  LITTLE  BIT  LARGER  THAN  EXPECTED  TO  PRPv-ENT 

13  REM  INSTABILITIES  AND  LIMITRIDING. 

14  REM 

15  REM  IN  SHORT  INSERT  VALUES  FOR 

16  REM  -      F(0)  AS  GKl),  LOWER  LIMIT  OF  F"  (0>  AS  G6  AhO 

17  REM  -      UPPER  LIMIT  OF  FMO)  AS  G7  IN  LINE  120 

IS   REM  -      LOWER  LIMIT  OF  NGI  AS  HG  AND  UPPER  LIMIT  OF 

19  REM  -      NGI  AS  H7  IN  LINE  130 

20  REM  -      WEIDGE  ANGLE  B  AND  FLUID  PROPERTIES  CONSTANT  R 

21  REM  -      IN  LINE  lAO 

22  REM 

23  REM  THE  OUTPUT  WILL  CONSIST  OF  TWO  TABLES  GIVING 

24  REM  FG,FG',FG",FL,FL'  AND  FL."  AS  A  FUNCTION  OF  THE 

25  REM  DIMENSIONLES  PARAMETERS  NG    AMD  NL, PLUS  A  HEADING 

26  REM  OF  THE  MOST  IMPORTANT  PARAMETERS  AND  CONSTANTS. 

27  REM 

28  REM 

100  DIM  G0( 101),  Gl(lOl),  G2<101),  G3(101> 

105  DIM  L0<  101),  LKlOl),  L2(101),  L3(101) 

110G0<1)=0:  Gl(l)=-.5:  Ga(l)=0:  G6=0:  G7=2 

115  H6=.95:  H7=1.05 

120  G4=G&:  G5=G7 

125  A=l:  B=l:  R=200 

130  Z=0:  Z1=0 

135  Ha=H6:  H9=H7 

140  H=(H8+H'9)/2 

145  GS=G6:  G9=G7 

150  G3<  l)  =  <GS-t-G3)/2 

155  IF  H=0  THEN  400 

160  IF  Z  <    35  THEN  ISO 
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165  PRINT  "AFTER  35  ITERATIONS  IN  GAS  UNABLE  TO  MEET  REG'S" 

170  PRINT  H,G1( 101),Ga< 101),G3< 1) 

175  GOTO  495 

ISO  H0=H/100 

185  FOR  1=1  TO  100 

190      H1=G2(I)*H0 

195      H2=G3(I)*H0 

200      H3=H0*(B*G2(  I>*G2(  I>-B-A-*C1(  I>*G3(  I  )  ) 

£05      G1<I  +  1)=G1<I)^-H1 

210      Ga(I+l)=G2(I)+H2 

215      G3<I+1)=G3<I)^H3 

220      G0<I+l)=G0<I)+HO 

225  NEXT  I 

230  Z=Z>1 

235  IF  ABS<G1(101))  <  .OOOOOOOl  THEN  255 

240      IF  Gl(lOl)  >  0  THEN  250 

245  Ga=G3(l>:  GOTO  150 

250  G9=G3(1>:  GOTO  150 

255L0<1)=0:  L1(1>=G1(101):L2(1)=G2(101>:  L3( 1 ) =G3( 101) /R 

260  M=10 

265  M0=M/100 

270  FOR  1=1  TO  100 

275      M1=L2(I>*M0 

280      M2=L3(I>*M0 

285      M3=M0^<B*t.2(I)*L2(I)-B-A*Ll(  I)*L3(I>  ) 

290      L1(I  +  1)=L1<I)^-M1 

295      L2<I-H)=L2(I)^'M2 

300      L3(I-H)=L3(  I)+M3 

305      L0<  I-H)=L0<  I  )+M0 

310  r-JEXT  I 

315  Z1=Z1+1 

320  IF  Zl<  35  THEN  335 

325  PRINT  "AFTER  35  ITERATIO^iS  IN  LIQ  UNABLE  TO  MEET  REG'S" 

330  GOTO  415 

335  C5=L2(101)-1 

340  IF  H=0  TI-EN  370 

345  IF  ABS<L3(101))  >  .OOOl  THEN  355 

350      IF  ABS<C5>  <  .OOOOl  THEN  410 

355  IF  C5  >  0  THEN  365 

360  H9=H:  G6=G3(1):  Z=0:  goto  140 

365  H8=H:  G7=G3(1>:  Z=0:  goto  140 

370  IF  AeS(L3(101))  >  .OOOl  THEN  380 

375      IF  ABS<C5)  <  .00001  THEN  410 

380  IF  C5  >  0  THEN  3'K) 

385  Ga=G3(l):  GOTO  395 

390  G9=G3(1):  GOTO  395 

3-35  G3<  1)  =  (G8-hG9)/2 

400  L3<  n=G3(l):  L0(1)=0:  L1<1)=G1(1):  L2tl)=L2(l> 

405  GOTO  260 

410  PRINT  "LAMINAR  FLOW  OVER  WEDGE  WITH  BLOWING" 

415  PRINT  :  PRINT  :  PRINT 
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450  PRINT  "8  =  ",B 

455  PRINT  "R  =  "  ,  R 

430  PRINT  "FG(0>       =  ", Gl< 1 )," INITIAL  LIMITS 

435  PRINT  '•FGDPRIME(0)=  "  ,  G3<  1  )  ,  G4,  G5 

440  PRINT  "NI  =  "  ,G0<101  ),H€.,H7 

445  PRINT  :  PRINT  :  PRINT 

450  PRINT  -NG", "FG", "FGPRIME", -FGDPRIME":  PRINT 

455  FOR  1=1  TO  101  STEP  10 

4€.0      PRINT  G0<I),Gl(I),Ga<I),G3<I) 

465  NEXT  I 

470  PRINT  :  PRINT  :  PRINT 

475  PRINT  "NL", "FL", "FLPRIME", "FLDPRIME":  PRINT 

480  FOR  1=1  TO  101  STEP  10 

485  PRINT   L0( I),L1(I),L2( I),L3< I) 

490   fvHXT    I 

495  PRINT  :  PRINT  :  PRINT  :  PRINT  "  "  "  "  "EhC" 
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APPENDIX  B 
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Figure  6.5   Dimensionless  velocity  profiles 
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Table  5.2   Comparison  Between  Single  Fluid  Flow 
and  Nickel  (abb  8)  [5] 


_s 

fq(0) 

'," 

(0) 

Program 

Nickel 

0 

0 

0.4365 

0.46 

0.25 

-.2 

0.5888 

0.60 

0.25 

-.5 

0.4399 

0.44 

0.25 

-1 

0.2644 

0.27 

0.5 

+1 

1.5758 

1.62 
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Table  6 . 1  Comparison  Between  V7all  Shear  Stress 
for  Large  Blowing  Parameter  and  Wall 
Shear  Stress  for  Single  Fluid 


6 

fq(0) 

f  '  • 

g 

(0) 

Program 

Nickel  [5] 

0 

-1 

0.0628 

0 

0.2 

-1 

0.2791 

0.25 

0.4 

-1 

0.4308 

0.40 

0.6 

-1 

0.5590 

0.53 

0.8 

-1 

0.6729 

0.66 

1.0 

-1 

0.7767 

0.76 

0.2 

-1.25 

0.1877 

0.18 

0.4 

-1.25 

0.3336 

0.33 

0.6 

-1.25 

0.4593 

0.45 

0.8 

-1.25 

0.5716 

0.57 
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Table  6 . 2   Local  Shear  Stress  Reduction  (t/t.^) 
as  a  Function  of  Wedge  Angle  and 
Blowing  Parameter  for  R  =  200 


f  (0) 

B 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

0.005 

0.01 

0.015 

0.02 

0.03 

0.04 

0.05 

0.075 

0.1 

0.15 

0.2 

0.35 

0.5 

0.75 

1.0 

0.508 

0.348 

0.265 

0.215 

0.155 

0.122 

0.100 

0.0691 

0.0526 

0.0353 

0-0262 

0.0138 

0.00827 

0.00337 

0.000720 

0.385 

0.252 

0.188 

0.150 

0.107 

0.0833 

0.0681 

0.0468 

0.0355 

0.0239 

0.0179 

0.00978 

0.00634 

0.00351 

0.00212 

0.331 

0.212 

0.156 

0.124 

0.0874 

0.0676 

0.0551 

0.0377 

0.0287 

0.0193 

0.0145 

0.00826 

0.00567 

0.00363 

0.00261 

0.299 

0.187 

0.137 

0.108 

0.0759 

0.0585 

0.0476 

0.0325 

0.0247 

0.0167 

0.0127 

0.00744 

0.00534 

0.00372 

0.00289 

0.276 

0.171 

0.124 

0.0972 

0.0681 

0.0524 

0.0426 

0.0291 

0.0221 

0.0150 

0.0115 

0.00692 

0.00515 

0.00379 

0.00309 

0.259 

0.158 

0.114 

0.0894 

0.0624 

0.0479 

0.0389 

0.0266 

0.0202 

0.0138 

0.0106 

0.00656 

0.00502 

0.00386 

0.00324 
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